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Abstract 



Two important notions of integrability for discrete mappings are algebraic integrability and 
singularity confinement, have been used for discrete mappings. Algebraic integrability is related to 
the existence of sufficiently many conserved quantities whereas singularity confinement is associated 
with the local analysis of singularities. In this paper, the relationship between these two notions is 
explored for birational autonomous mappings. Two types of results are obtained: first, algebraically 
integrable mappings are shown to have the singularity confinement property. Second, a proof of the 
non-existence of algebraic conserved quantities of discrete systems based on the lack of confinement 
property is given. 

I INTRODUCTION 

One of the first notion of integrability was introduced by Liouville in the nineteenth century for 
Hamiltonian systems in classical mechanics. This type of integrability is based on quantities whose 
value do not change in time, the so-called constants of the motion or first integrals. For instance, 
in many mechanical systems the total energy and the linear and angular momenta are conserved 
through the dynamics and are examples of such constants. Liouville's fundamental contribution was 
to prove that if a given Hamiltonian systems admits enough first integrals, it can be solved explicitly 
by quadratures. More precisely, a 2n-dimensional Hamiltonian system is said to be Liouville-integrable 
if it admits n functionally independent first integrals in involution (that is their Poisson Brackets 
commute) and such systems can be integrated by quadrature and enjoy a particularly simple topology 
(the flow lives on products of tori and cylinders). For systems of n first-order ordinary differential 
equations (ODEs) algebraic integrability is defined as the existence (n — 1) functionally independent 
first integrals that are algebraic functions of the dependent variables. The existence of such first 
integrals is important for integrability as they can be used to reduce the dimensionality of the system. 

In the particular case of Hamiltonian systems, algebraic integrability coincides with the notion 
of superintegrability. Indeed, for a 2n-dimensional Hamiltonian system, the 2n — I first integrals 
needed for the property of algebraic integrability are much more than the n ones needed for complete 
Liouville- Arnold integrability. 

A different notion of integrability was introduced by Painleve in the beginning of the 20th century 
|21[I5]. Although the aim of Painleve was not to define a notion of integrability but rather to build 
new functions, his property is today widely used for the detection of integrable systems. An ODE 
is said to possess the Painleve property pi ll7l l()j if its general solution is single-valued in its maximal 
domain of analytic continuation. The restrictions the Painleve property impose on the solutions are 
so strong that an ODE exhibiting it may be considered for all practical purpose integrable. Many 
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formal links between the Painleve property and other notions of integrability have been established 
(see for example 013 )■ Despite the fact that there is no general algorithmic way to obtain sufficient 
conditions for a given ODE to have the Painleve property, necessary conditions can be derived. The 
algorithmic procedure to obtain such conditions is known as the Painleve test or singularity analysis 
and is based on a local analysis of the solutions around movable isolated singularities. 

The notions of algebraic integrability and singularity analysis of the solutions can be extended to 
finite difference equations. One of the most effective way to perform singularity analysis in this context 
is given by singularity confinement "T*'!! , T2",'53"'T,'25"26', 13 , 14 , f 5",T6' which, like the Painleve test for 
ODEs, imposes conditions on the singularities of the solutions. Although there exist many striking 
particular examples which indicate that singularity confinement is closely related to other notions of 
integrability, there is no clear-cut result which establishes this relationship in a formal and general 
framework. This paper explores the relationship between the notion of algebraic integrability and the 
property of singularity confinement in discrete mappings. 

In the case of systems of ODEs, Yoshida |27ll28llH] proved that the degree of a rational first integral 
is closely related to some exponents that can be obtained through the application of the Painleve 
test. In this paper, we establish an equivalent result for discrete systems by showing that singularity 
confinement analysis provides necessary conditions for the existence of an algebraic first integral. 
Moreover, a lower bound on the degree of rational first integrals can be obtained. 

The rest of this article is organized as follows. In Section m we define the notion of singularity 
and singularity confinement for birational autonomous discrete dynamical systems. In Section IIIII 
we study the local consequences of the existence of algebraic first integral on the confinement of the 
singularities in the two-dimensional case. In Section IIVI we extend the results of Sections IIIII to 
arbitrary dimensions. Applications are discussed in Section Ivl 



II FORMULATION OF THE PROBLEM 

In this section we introduce and illustrate the notion of singularity confinement by considering a 
particular class of two-dimensional autonomous dynamical systems. These examples are used to 
motivate a formal definition valid in arbitrary dimensions. 

Throughout the article the following notation will be used: If g is an analytic function from C 
to C, then Dg(x) is its Jacobian matrix evaluated at x G C^. When p = r, the Jacobian is denoted 
by det(Dg(x)). If / is a complex-valued rational function on C, then num(/) and den(/) denote 
respectively the numerator and denominator of /. 



II. 1 Simple examples 

The basic idea of singularity confinement is to consider the properties of solutions close to some 
singularities. To illustrate this concept, consider the class of two-dimensional complex dynamical 
systems of the form 

7+0=f(.„,,„).(^(^-^"M, (2.1) 
where 5 is a complex- valued rational function on C^, 



q{x,y) 



and p, q are relatively prime polynomials. Two types of singularities can be distinguished for such 
systems. The singularities of first type are the roots of q in C^. That is the values (a;, y) at which the 
vector field t{x,y) in l|2.1() has a singularity. The singularities of second type are the points where the 
Jacobian of f is zero (this type of singularities was first considered in |12|). As an example, consider 
the following discrete dynamical system: 



^ =f(a;„,y„)= ( ^" ^" + + a!l | , (2.2) 



Vn+l 
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where a and b are complex numbers. The analysis is performed for generic values of the parameters a 
and b. Clearly the right-hand-side of (|2.2(l presents a singularity if Xn = 0. To study this singularity 
we study the behavior of nearby solutions by introducing a small complex perturbation e, with |e| > 0. 
We then set ya to be any complex number and xq ^ e. The iterates Xp and yp are then determined for 
any positive nonzero integer p by the discrete dynamical system (|2.2() . After one iteration, we have 

xi = -^+a-yo + 0{e). (2.3) 
In the limit e — > 0, jxij — > cx) and, using H2.2|l again, 

X2 = -^+ya + 0{e), (2.4) 

X3 = -e + 0(e^). (2.5) 

The limit e ^ of the next iterate X4 is well defined and given by yQ. Despite the fact that the function 
f defined by (|2.2|l is not well defined on any point of of the form {Q,y), the limit {x,y) (0, y) of 
{'^{x,y) exists and is given by 

lim,„ {\x,y)^( y). (2.6) 



{x.,v)~.{0.y) ' V 

Moreover, since the value of the limit (|2.6(l depends on y, the same limit applied on the Jacobian of 
f* is nonzero. The singularity is thus "confined" between the iterates and 4 and does not propagate 
further. This property of confinement is very particular and we show in this paper that it is closely 
related the existence of a first integral for the system. A first integral for a discrete dynamical system 
of the form (|2.1|) is defined to be an analytic non-constant complex-valued function / defined almost 
everywhere on which is preserved by f , that is 

/(f(x))=/(x), (2.7) 

for every x e for which H2.7() makes sense. It is important to note that if f is not defined at a point 
X* where / is well-defined, then the value of / at x* is still preserved under f in the limit sense, i.e. 

lim /(f(x)) =/(x*). (2.8) 

X — >x* 

This is a consequence of the fact that f and / are continuous and defined almost everywhere in C^. In 
the particular case of system ()2.2|l , one can check that the following polynomial 

I{x,y) = xy{x + y - a) - b{x + y) (2.9) 

is a first integral. Moreover, following (|2.8|) . the first integral / is preserved in the limit sense through 
the iterations of the singularity of (|2.2|l . That is, 

-6y = 7(0,2/) = lim I o {r{x,y)) , t = 1,2,3,4. (2.10) 

Note that most integrable systems used in this paper are particular cases of the so-called QRT-family 
of mappings P^IT^ . 

The confinement property for the solutions of (|2.1() does not hold in general. Although there is no 
formal result establishing how rare this confinement property might be, the fact that it does not hold 
for most systems of the form H2.1|l has been well established in the literature jllll23ll^l2()lll3lll4j . 

The following system is an example where singularities are not confined: 

b 

) = f{Xn,yn) =1 Xl \ ■ (2.11) 

yn+l 



Consider the singularity {x*,y*) = {0,y) and, as before, introduce xq — e and yo C nonvanishing. 
The next iterates are xi ^ b/e^ + {a — yo) — e, X2 = —bje' + yo + 0{e'), and 2:3 = e + ©(e^). Further, 
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we find that, unlike the previous example, the appropriate cancellations allowing for confinement do 
not occur and the next iterates are also diverging at e = 0, a;4 = 2b/e^ + (a — yo) — e + 0{e^), and 
X5 = —2b/e^ + yo + C'(e^). In general, the sequence of limits e ^ of the |a;p| for p > repeats the 
formal triplet (00, 00, 0) indefinitely. This assertion is proven by induction based on the equality 

-{n+l)b/e^ 




{n + l)b/e^ + (a 




(2.12) 



for any positive integer n. Therefore, the singularity is not confined. It is generally believed that 
systems which lack the confinement property will not be integrable. However, there is no definite 
result attached to this belief. We prove in the next sections that we can use the analysis of the 
dynamics nearby the singularities can be used to conclude that H2.11|l does not admit an algebraic first 
integral. 

The following example illustrates the occurrence of the second type of singularity. Consider 

^ -nx„,y„)= i + I . (2.13) 



Vn+l 



Xn 



At first sight, the point (x„,0) seems to be the only singularity. However, the Jacobian of f in H2.13|l 
vanishes whenever a;„ = 1. If a;o 1 and yo ^ 0, the next iterate as determined by H2.13|l is given 
by xi ~ 2. We then have X2 = 5, X3 — 12. It is not difficult to see that Xi grows with respect to i 
and that all iterates, of (1, y) under H2.13|l with y ^ are independent of y. The singularity here only 
appears in the inverse of the Jacobian matrix and is not confined. 

In the light of the previous examples, we can define the confinement property for singularities of 
discrete dynamical systems of the form ()2.1|l . A general definition will be given in the next section. 
A singularity of the dynamical system (|2.1|) is defined to be any point {x*,y*) in C'^ at which the 
right-hand-side of H2.1() is undefined or at which the Jacobian of f is zero. A singularity (x* , y*) is said 
to be confined if there exists a positive integer N such that both limits 

lim f^(x,?/), lim det(Df^(a;,y)) (2.14) 

exist and the second limit is nonzero. The smallest number N having this property is referred to as 
the confinement number, that is the number of steps necessary for confinement. 

The next example illustrates the case when a singularity of the second type is confined. 

/X / {xn - l/a)ixn ~ a) \ 

( ^^+^1 j = iixn,yn) = ynixn ~ 6) (x„ - 1/6) , (2.15) 

\ Xn J 

where a and b are two nonzero distinct complex numbers which are also both distinct from 1/a and 
1/6. The right-hand side (RHS) of (|2.15|) presents a singularity of the second type if x„ — a because 
the Jacobian of f vanishes, li xq — a and t/o 7^ 0, the next iterate as determined by (|2.15() will be 
xi = 0, then X2 = ^/a and x^ takes an indeterminate form 0/0. Again, we introduce a small complex 
perturbation e, |e| > 0, then set xq = e to find that when e — > 0, 

1 

x^^yo + a , (2.16) 

a 

and we conclude that the singularity is confined. The dynamical system (I2.15|l has another singularity 
of the first type if x„ — 6. Applying the same procedure as before with a perturbation e, one finds 
that the limit e — > of |a;i|, X2 and x^ gives, respectively, oo, 1/6 and byo/{b + yo{l — 6^)). Thus, any 
singularity of the form {b,yo) with yo ^ 6/(6^ — 1) is confined in 3 steps. System H2.15|l admits the 
following rational first integral 



Singularity Confinement and Algebraic Integrability 



5 



The following example shows that the existence of a first integral does not imply that all singularities 
are confined. The system 

^ =f(x„,y„)= f V (2.18) 



has the following first integral 

Iix,y)^^^. (2.19) 

The RHS of (|2.18|) has a pole of order 1 at any point of of the form (x, 0) with x 0. As before, if 
we study this singularity by perturbation we conclude that there is no confinement since the Laurent 
expansions of the iterations of {x, e) under the discrete system H2.18|l around e = are all divergent. 
Although this singularity is not confined, one notices that any point of the form {x, 0) with x ^ does 
not have a unique preimage under H2.18|l . Indeed, the expression defining given by 



{-Hx,y)=\ (2.20) 




has a singularity of second type at (a;, 0) and all the points of this form are mapped to (0, 0) under f 
This example shows that if one wishes to establish a relationship between singularity confinement and 
first integrals, particular care should be given in distinguishing different singularities. 

II. 2 Preliminaries 

The notion of singularity confinement will now be defined in a more general setting. We consider 
p-dimensional autonomous discrete dynamical systems described by birational mappings on C^. A 
birational mapping on is a C-valued rational function which is one-to-one almost everywhere in 
and whose inverse, where it exists, is also represented by a rational function. The corresponding 
discrete dynamical systems associated are written: 

x„+i = f(x„), x„ e C^ n e Z, (2.21) 

where f is a birational mapping on with no explicit dependence on n. The ith iterate of xq G C 
under H2.21(l is said to exist if the following limit 

lim r(x) (2.22) 

x^xo 

exists. The ith iterate is then defined to be that limit and is denoted F*(xo). 

Definition. A singularity of type I for the discrete dynamical system x„+i = f(x„) is a point 
x* G at which f is not defined. The set of singularity of type I associated with the mapping f in C 
is denoted Si{{). 

Definition. A singularity of type II for the discrete dynamical system Hn+i — f(x„) is a point 
x* e such that det (Df (x*)) = 0. The set of singularity of type II in C is denoted S'//(f). 

In the particular case of birational mappings, we have 

Si{f) = |x* e CP dcn(/,)(x*) = 0, for some i G {1,2, . . . , (2.23) 

Sii{i) = {x* GCf det(Df(x*)) = o}. (2.24) 

We call the sets ^/(f) and Sn{f), the singular sets of first and second type, respectively. The singular 
set S{f) is defined as the union of both sets of first and second type 

S{i) = Si{i)USii{i). (2.25) 
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For both types of singularity we can define the property of confinement. 

Definition. Let x* G C 6e a singularity of type I or II for the system x„+i = f(x„). The singularity 
is said to be confined if, for some ic G N^", the iterate F''^(x*) exists and linix^x* dot (_Df''^(x)) ^ 0. 

The lowest such ic is referred to as the confinement number that is the number of steps necessary 
for confinement. From now on, we assume that f is a birational mapping. The following lemma proves 
to be useful. 

Lemma 2.1. The image o/5//(f) under the mapping f lies inside an algebraic variety of codimension 
2. 

That is, the birational mapping f is one-to-one only on the subset of defined by 

Of = CP\5(f). (2.26) 
As a consequence, the function f restricted to Of defines a rational diffeomorphism 

f : Of ^ O^ = f (Of), (2.27) 
whose inverse is also a rational function. 

Proof. Assume that ^//(f) ^ 0. Let xq G Su{^) and let xi = f(xo). From the inverse function 
theorem, it follows that 

lim II det(Df"i(x)) II =oo. (2.28) 

Thus, for some j's in {1, 2, . . . the denominator of (f~^)i evaluated at Xi is zero. Denote those i 
by ifc. A: = 1, 2, ... , r. We now prove that at least for one ifc, num also vanishes at Xi. By 

contradiction, suppose that none of these numerators vanishes at xi . Thus 

lim ||f"^(x)|| = oo. (2.29) 

On the other hand, we have 

hm ||f-i(f(x))|| =xo, (2.30) 

X-+XO 

which contradicts equation (|2.29|) . Thus, xi must be in one of the following sets 

A, = {x e C^l num((f-i),)(x) = den((f-i),)(x) =0}, le {f , 2, . . . (2.31) 

This concludes the proof since each set Ai is defined as the intersection of the zero level sets of two 
relatively prime polynomials. □ 

As mentioned earlier, the singularity confinement turns out to be closely related to the existence 
of a first integral. This relationship will be explored in the next section. We define the notion of a 
first integral. 

Definition. A first integral o/ x„+i = f(x„) is an analytic function I : U C C ^ C where U 
is a dense subset of such that 



(a) (dL^OL ai_ 

* ^ V OXl ' OX2 ' ' OXp 

(b) /(f(x))=/(x). 



7^ almost everywhere in U . 



A first integral can exist in several forms. 

Definition. A first integral I o/ x„+i — f(x„) is polynomial (resp. rational^ if the function 
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/(x) is a polynomial (resp. rational) function ofx G C^. A complex-valued function /(x) is algebraic 
over C if there exist s > and qq, . . . ,qs rational in x such that 

go + gi/ + --- + 9.r-0. (2.32) 

// s is the smallest positive integer such that 1^2.!^!^} holds, the relation Y2.3^) is referred to as the 
minimal polynomial of f . A first integral I is algebraic if the function /(x) is an algebraic function. 

Many notions of integrability are used in the literature. In this paper, we concentrate on algebraic 
integrability. In order to introduce the latter, we are interested in determining precisely "how many" 
algebraic first integrals a given discrete system must admit to be algebraically integrable. To further 
clarify this problem, suppose that the algebraic functions (x), v42(x), . . . ,Ar{x) defined on are 
all first integrals for the system (|2.21|) and F{zi, Z2, . . . , Zr) is an algebraic function on C". Then, it is 
clear that A(x) = F{Ai{x), A2(x), . . . , Ar{x.)) is also a first integral. However, this new first integral 
adds no knowledge to the given problem since it depends on the other first integrals. 

Definition. Let yli (x), yl2(x), . . . , A,.(x) be smooth complex-valued functions defined on a domain 
D dCP. Then 

(a) Ai, A2, . . . ,Ar are functionally dependent if, for each x G D, there is a neighborhood U of x 
and a smooth complex-valued function F{zi, Z2, ■ ■ ■ , Zr) not identically zero on any subset of C 
such that 

F(Ai(x), A2(x), . . . , A(x)) = 0. (2.33) 

(b) Ai, A2, ■ ■ ■ , Ar are functionally independent if they are not functionally dependent when 
restricted to any open subset of D. 

A simple way to determine if a set of functions are functionally independent is given by the follow- 
ing theorem (see for example 20 ). 

Theorem. Let A = {Ai, A2, . . . , A^) be a set of smooth functions from a domain D ^ VP to 
C. Then, Ai(x), A2(x), . . . , Ar(x) are functionally dependent if and only if DA(x*) has rank strictly 
less than r for all x* G _D. 

We introduce the notion of integrability considered in this paper. 

Definition. The discrete dynamical system x„+i — f(x„) is said to be algebraically integrable 

if it admits {p — 1) algebraic functionally independent first integrals. 

The following theorem is a direct extension of Brun's [5l ll9) theorem to discrete systems and is useful 
to extend the results we obtain for rational first integrals to algebraic ones. 

Theorem.: // a discrete dynamical system o/x„+i f(x„) has k functionally independent algebraic 
first integrals, then it has k functionally independent rational first integrals. 

Proof. Let / be an algebraic first integral and let 

P{I) ^ go + qi/ + . . . + + I' (2.34) 

be its minimal polynomial, where qi is a rational function of x. Since / depends non trivially on x, let 
i be such that is a non constant rational function of x. Since /(f(x)) = -/^(x), we have that 

go(f(x)) - go(x) + (<7i(f(x)) - -7i(x))/(x) + • • ■ + (g.-i(f(x)) - g,_i(x))r-i(x) = 0. (2.35) 

Since P is minimal, we have (7i(f(x))) = (J'j(x) and each qi is a first integral. Now, let /' be another 
independent first integral whose minimal polynomial is 

P(/') =q'^ + q[l' + ■■■ + q:,_,r'-' + l'^' . (2.36) 
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The independence between the two first integrals imphes that there exist i < s, i' < s' such that qi, qii 
are two independent non constant rational first integrals. By induction, one can build k independent 
rational first integrals. □ 

Therefore, rational integrability implies algebraic integrability and it is sufficient to consider discrete 
dynamical systems of the form 12.21|l which admit (p — 1) functionally independent rational first 
integrals Ri (x) = P, (x) / Qi (x) where i = l,2,...,p— 1 and Pi and Qi are relatively prime polynomials . 
Let Li^c be the level set of the first integral Ri corresponding to the value Ci G C. Since f is continuous, 
it leaves invariant the closure of the level set given by 

r,,c = {x e CP |P,(x) - cQ,(x) = 0}. (2.37) 



III The two-dimensional case 

In this section we study the local implications of algebraic integrability on the confinement of 
singularities for two-dimensional discrete dynamical systems. We consider a discrete dynamical systems 
of the form H2.21|l which is algebraically integrable, that is, it possesses a non trivial rational first 
integral i?(x) = P(x)/(5(x) where P and Q are some polynomials over C^. Let Lc be the level set of 
R corresponding to the value c G C. Then, as noted before, f leaves invariant the closure of the level 
set Lc- 

We first focus our attention on singularities of the first type. The singular set of first type S'/(f) is 
an algebraic variety since it is defined as the set of zeros of polynomials. We consider its irreducible 
decomposition 

d 

Si{i)^[jsf{i) (3.1) 

i=l 

where each irreducible component s'f\i) is defined as the zero set of an irreducible polynomial on C^. 

In Example (|2.18|1 . we showed that the existence of a rational first integral does not ensure that all 
singularities are confined. However, in this case it is not possible to "enter" the singularity, meaning 
that it does not have a unique preimage under the mapping. Moreover, we know that confinement 
is a generic property, that is, it is only concerned with dense subsets of irreducible components of 
Sj. For instance, in Example (|2.15(l all singularities of type I of the form {b,yo) confine in 3 steps 
except a particular singularity defined by yo = b/{b'^ — 1). We state the general theorem relating the 
confinement of singularities to the existence of a rational first integral. 

Theorem 3.1. Consider a two-dimensional birational mapping x„_|_i = f(x„) and assume it has 
a rational first integral P(x) = P(x)/(5(x). Then, for each irreducible component Sj''\{), we have 
either 

(a) There exists kj g N^*^ such that 

lim f'=^(x) (3.2) 

X — ■ X " 

exists for almost all x* G Sj''\i), or, 

(b) There exists Ij G N^*^ such that Sj''\{) lies inside the singular set of the birational mapping 
defining f . 

Proof. Without loss of generality, take j = 1. Consider two cases: 
1) S^/\i) the forallcGC. 

Let X* G 5|^''(f) such that limx^x* ||/j(^)ll exists or is infinity (as opposed to undefined) for 
i = 1,2. Suppose that x* G Lc* for only one c* G C and that there is a neighborhood of x* 
on Lc* whose only intersection with any of the singular sets associated with the mappings P, 
i = 1,2, . . . ,q+l (q is defined below) is x* itself. The set of points satisfying the above properties 
is dense in S'j^^(f). Let q be the number of paths Pi(e), i — 1,2, ..,q in such that 
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(i) Pi(e) G ic* for e small enough. That is, 

P(p,(e))-c*g(p,(e)) = 0. (3.3) 

(ii) ||pi(e)|| — > oo as e ^ and at least one of the components of Pi(e) is of the form 1/e. 

Note that the number q can depend on the complex number c* but the values it can take for 
different c* can be bounded above by a number depending only on the degrees of P and Q. 
Now, consider another path x(ct) G Lc* Vu such that 

limx(cr)=x*. (3.4) 

cr— >0 

Since there is a neighborhood of x* on Lc* whose only intersection with S{i) is x*, f(x(cr)) is 
well-defined for a nonzero and small enough. Moreover, since f preserves Lc*i f(c) also lies in 

Lc' , and 

lim||f(x(a))||=^. (3.5) 

Therefore, there exists a change of variable a = <T{e) of order 0(6*) for some s > 0, so that 
f (x(cr)) is one of the paths pi defined above: 

f(x(a(e)))=p,(e) (3.6) 

for some 1 < i < q. Our claim is that there exists m with m < q + 1 such that the limit as 
cr of f'"(x((T)) exists and is finite. Suppose by contradiction that the iterates f*''(x(CT)) 
k = 1, 2, 3, . . . , g + 1 are all divergent at cr = 0. Each iterate f'^(x(cr)) is associated with a unique 
path pij. (e). This implies that the k^^ preimage of p^^ (e) is finite in the limit e — > and no other 
f'(x((T)) for l<k<l<q+l corresponds to the same pi^{e). Hence the relation between the 
q+1 iterates that are divergent at cr = and the paths Pi(e) is injective. This is a contradiction 
since there are only q paths Pi(e). 

There exists a c G C such that S'|^''(f) C Lc- Let 



L,= U4^^ (3.7) 



be the irreducible decomposition of Lc for any c G C. If 5|^''(f) lies inside the singular set of 
, then the theorem is proved. Otherwise, let k G N>° be such that S'j^^ (f ) does not he m a 
singular set associated with the mappings f ~™, m — 1,2, . . . ,k. Then, 



f-m(^(l)(f)) ^ U l<l<k (3.8) 



where each B„i is a subset of {1,2, ...,r}. The k sets in (|3.8|l are distinct, otherwise, by 
contradiction there exist mi,TO2 G such that 



f-'"i(S'j^'(f)) =f-™2(S'j^^(f)), l<mi<m2<fc. (3.9) 



This implies that 



S'W(f) = f-i (f(™i-™2+i)(5(i)(f))^ (3.10) 

which is a contradiction since f is undefined almost everywhere on 5|^^(f). Hence the k sets in 
(|3.8|l are distinct. However, there are only finitely many possibilities for the sets Bm in H3.8|l . 
Hence S'|^''(f) must lie inside the singular set of f^'^ for some h. 

□ 
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Notice that the condition H3.2|l does not mean that the singularity is confined since the condition 
that the Jacobian be nonzero must be satisfied. In order to solve this problem, the singularities of 
type II have to be considered. The singular set of second type is not, in general, an algebraic variety 
but its closure 

SVRf) = {x* e C num (det (Df (x*))) = 0} (3.11) 
is. Hence, it is possible to introduce the irreducible decomposition 



Su{i) = \JS^i>{{), (3.12) 

where each s\j\t) is a subset of S'//(f). It is now possible to state the general theorem for singularities 
of second type. 

Theorem 3.2. Consider a two-dimensional birational mapping x„-|_i = f(x„) and assume it has 

a rational first integral R{x) = P(x)/Q(x). Then, for each irreducible component s\j\{), we have 
either 

(a) There exists kj G such that 

lim det (Df'^J (x)) (3.13) 

exists and is nonzero for almost all x G Sjp{{), or, 

(b) There exists Ij £ such that Sjp{{) lies inside the singular set o/f^'^ . 

Proof. The proof is similar to the proof of Theorem 3.1. Again, without loss of generality, we take 
j = 1 and consider two cases. 



1) S^jj\i) (^Lc foraUceC. 

The image ofS\j\f) under the mapping f is a finite subset of C^, thanks to Lemma 2.1. Moreover, 
Sjj\f) is not contained inside any level set. For any point in the image of Sjj\{) it is therefore 
possible to choose two distinct preimages in s\j\f) which lie in two different level sets. Since f 
preserves Lc, every point in the image of Sjj\i) must lie in the set 

P|L, = {xeC2|p(x)=Q(x)=0}. (3.14) 

cGC 

The elements of this finite set are denoted Si,i — 1,2, ... ,m. Note that in the case of a polynomial 
first integral, Q(x) = 1 and this set is empty. 

Now, choose x* G s[j\f) such that x* S Lc* for only one c* G C and there is a neighborhood 
of X* on Lc' whose only intersection with any of the singular sets associated with the mappings 
f*, i — l,2,...,g + r + l {q and r are defined below) is x* itself. The set of points satisfying the 
above properties is dense in S^j^\{). 

Consider the q paths Pi(e) described in the proof of Theorem 3.1. Consider also the r paths 
Wi(e), i = 1,2, ..,r in satisfying 

(i) Wi(e) e Lc' for e small enough. 

(ii) ||wi(e)|| Sj. as e ^ and one of the components (fc = 1 or 2) of Wi(e) is of the form 
(sj. )fc + e for some 1 < ji < m. 

Now, consider another path x(cr) G Lc* such that 

limx(cr)=x*. (3.15) 

Since there is a neighborhood of x* on Lc* whose only intersection with S{i) is x*, f(x(cr)) is 
well-defined for a nonzero and small enough. Moreover since f preserves Lc*, f(x(a)) G Lc*, and 



lim f(x((T)) = Si 



(3.16) 
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for some i such that 1 < i < r. FoUowing the proof of Theorem 3.1, f (x((t)) can be associated 
with a path Wj(e) (whose hmit as a ^ is s^). Suppose by contradiction that the d + r -\- 1 
iterates f'^(x(cr)) k = 1,2, ...,d + 7- + l are either divergent at cr = or have a hmit a ^ 
equal to one of the s^. Each iterate is associated with a path Pi^{e) or Wi^(e). This imphes 
that the unique k^^ preimage of p^^ (e) or Wi^, (e) must converge to x* as e — » 0. The preimage 
is unique because there is a neighborhood of x* on Lc» whose only intersection with any of the 
singular sets associated with the mappings P , i = 1,2, . . . ,q + r + 1 {q and r are defined below) 
is X* itself. Since x* is not one of the s^, there cannot be two iterates f'^^(x(cr)) and f'^^(x(cr)), 
ki, k2 < q + r + 1 associated with the same path. Since there are only q + r + 1 paths Pi(e) and 
Wi(e), there exists ak<q + r + l such that the limit cr — > of the fc"^ iterate of x(o') must 
be finite and not equal to one as cr ^ 0. Therefore, the limit cr —> of the Jacobian of f'^^ 
evaluated at x(cr) must be a finite nonvanishing number. This ends the proof in the first case. 

2) There exists a c G C such that Sj^\f) C Lc- The argument leading to the conclusion that s'fp {{) 
lies inside the singular set of f^'^ for some li is identical to the second case considered in the 
proof of the preceding theorem. 

□ 

With Theorem 3.2, it is now easy to see that singularities of type I satisfying the condition H3.2|l are 
confined. Indeed, suppose kj is the number considered in the first part of Theorem 3.1. If the Jacobian 

of f*^^ is zero almost everywhere on Sj(f), then s\''\i) becomes a subset of Sii{{''j ). According to the 
proof of Theorem 3.1 we know that S'j"'' (f ) ^ Lc for all c G C. Thus, as a consequence of Theorem 3.2, 
the singularities in s'f\{) must be confined. We can state a general theorem including both types of 
singularities. To do so, we introduce the irreducible decomposition associated with S'(f) 

d+s 

^(f) = IJ S'«(f) (3.17) 
where d and s are defined in H3.1|l and H3.12|) . 

Theorem 3.3. Consider a two-dimensional birational mapping x„+i = f(x„) and assume it has 
a rational first integral i?(x) = P(x)/Q(x). Let 5(f) be its singular set with irreducible components 
5'(j)(f)_ Then, for each j , either 

(a) Almost all singularities of S^-'\{) are confined in the same number of steps, or, 

(b) There exists Ij G such that ^'^'(f) lies inside the singular set of . 

IV Arbitrary dimensions 

In this section, we extend the result of the previous section to birational mappings of arbitrary 
dimension p. The mappings are assumed to be algebraically integrable, that is, they admit {p ~ I) 
functionally independent first integrals i?i(x) = Pi{x.)/Qi{x), i = 1,2, . . . ,p — 1. Then, the closure of 
the level sets associated with c G C^"^ is given by 

_ p-i_ 

Lc - {x G C^'lF,(x) - c, g,(x) =0, for 1 < i < p - 1} = fl L,,e,, (4.1) 

i=l 

with Li, a defined in (|2.37() . A set of constants c G is said to be regular for the first integrals 

R = (i?i, i?2, . . • , Rp-i) if DR(x) has rank equal to {p — 1) almost everywhere on R~^(c). The set Lc 
is then an algebraic variety of codimension {p — 1). 

Our main Theorems 3.1, 3.2, and 3.3 can be readily generalized to arbitrary dimensions and we 
only give here an outline of the proofs. 
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Theorem 4.1. Consider a p- dimensional birational mapping x„_|_i = f(x„) and assume it is al- 
gebraically integrable with first integrals i?i(x) = Pj(x)/(5j(x), i = 1,2, . . . ,p — 1. Then for each 

irreducible component sf ^t), we have either 

(a) There exists kj G such that 

lim f'=^(x) (4.2) 

X — 

exists for almost all x G ^^"'■'(f), or, 

(b) There exists Ij € such that Sj''\f) lies inside the singular set of the birational mapping 
defining f . 

Proof. Without loss of generahty, take j = 1. Consider two cases: 

1) S'j^\i) L^^c for all c G C^-^ and 1 < i < p - 1. 

Let X* G <S'}^^(f) such that linix^x* exists or is infinite (as opposed to undefined) for 

i = 1, . . . ,p. Suppose that x* G L^.* for only one c* G CP~^ which is regular and that there is 
a neighborhood of x* on Lc* whose only intersection with the singular sets associated with the 
mappings P , i = 1,2, . . . ,q + 1 {q is defined below) is x* itself. The set of points satisfying the 
above properties is dense in S'j^^(f). Let q be the number of paths Pi(e), i = 1,2, .., q in C^* such 
that 

(i) Pi(e) G Lc* for e small enough. That is, 

Pk{Pi{e))-clQk{Pi{e))=0, k = 1,2, . . . ,p - 1. (4.3) 

(ii) ||pi(e)|| ^ 00 as e — > and at least one of the components of Pi(e) is of the form 1/e. 

Using the fact that Lc» is of codimension p — 1, the rest of the proof follows exactly the similar 

case in the proof of Theorem 3.1. 

2) There exists a c G C and al<i<p— 1 such that Sj^\f) C -Li,c- 

This part of the proof is identical to the corresponding part in the proof of Theorem 3.1 except 
that Li^c is considered instead of Lc- 

□ 

Theorem 4.2. Consider a p- dimensional birational mapping x„+i = f(x„) and assume it is 
algebraically integrable with rational first integrals -Ri(x) = Pj(x)/Qj(x), i = 1, 2, . . . ,p — 1. Then, 

for each irreducible component Sjp{f), we have either 

(a) There exists kj G N^" such that 

lim det (i?f'^J(x)) (4.4) 

exists and is nonzero for almost all x G 5*1^ ''(f), or, 

(b) There exists Ij G such that s\p {{) lies inside the singular set o/f^'^. 

Proof. The proof is similar to the proof of Theorem 3.2 and the only difference comes in the first case, 
when Sjj\{) <^ Li^c for all c G C and 1 < i < p — 1. Following the proof of Theorem 3.2, for any 
point in the image of Sjj\f) it is possible to choose two distinct preimages in Sjj\f) which lie in two 
diflferent level sets. Indeed, from Lemma 2.1, the image of 5})^(f) under the mapping f lies in a subset 
of codimension 2 inside C^. Since f preserves ij_c, the points in the image of Sjp{{) must lie in the 
set 

fl Ii,e = {xGC2|Pi(x)=Qi(x)=0}. (4.5) 

ceCA<p-i 

Since the first integrals are functionally independent, the set defined above is finite. The rest of the 
proof is similar to the proof of Theorem 3.2 except that one should specify that c* has to be regular. □ 
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Theorem 4.3. Consider a p-dimensional birational mapping x„_|_i — f(x„) and assume it is 
algebraically integrable. Let S{i) be its singular set with irreducible components 5*^^ Then, for 
each j , either 

(a) Almost all singularities of S^^^{{) are confined in the same number of steps, or, 

(b) There exists Ij S N^'' such that S^^^{{) lies inside the singular set of i^^^ . 



V Corollaries and applications 

In the previous two sections, the existence of first integrals was assumed to find local information 
on the confinement property. Here, the information given by the singularity confinement property 
is used to obtain global information on the discrete dynamical systems such as the non-existence of 
algebraic first integral and the degree of possible rational first integrals. The results of this section 
can be considered as a discrete analog of Yoshida's Theorem |27II28[ I^ which, for ODEs, relates the 
Kovalevskaya exponents given by the Painleve test to the degree of a rational first integral. The first 
corollary is a direct consequence of Theorem 4.3. 

Corollary 5.1. Consider a p-dimensional birational mapping x„+i — f (x„) with a non-empty singular 
set S{{). If there exists an irreducible components S^^'^{f) in which almost all singularities are not con- 
fining and such that S^^^{i) (t S'(f~'^') for any positive k, then the system is not algebraically integrable. 

As an example of this corollary, we show that System 1)2.1 f |l . defined by 

f{x,y)=(^-'^-y + ^+^y^ (5.1) 

is not algebraically integrable. Recall that singularities of the form (0, y) are non confining. In order 
to use Corollary 5.1 and prove that this mapping does not admit an algebraic first integral, one must 
show that points of the form (0,2/) are not generically in the singular set of f~' for some I > 0. The 
singular set of f^^ consists of points of the form (a:,0). Thus, we must prove that generic iterates of 
(0, y) under f do not belong to the set of points of the form (x, 0). The inverse of the mapping H5.I|1 
is given by 

{-\x,y) = ( b ] . (5.2) 

\ -x-y + a+^ j ^ ^ 

This mapping can be obtained from f of (|5.I|) by interchanging x and y. Therefore, we can follow the 
analysis performed in H2.I2|I to perform a Laurent expansion in y to show that 



^ ^b/yl + xo + 0{yl) )• ^^■^> 



Moreover, we have 

nb/yi + {a- xo) - yo + C(2/o) \ / (^ + l)V2/o + (a - a^o) - Vo + 0(y§) 



F-3 



(5.4) 



-nb/y^ + xo + 0{y^) J \ -{n + l)b/y^ + xo + 0{y^„) 



where n is any positive integer. Hence, points of the form (0, y) are not generically sent to points of the 
form {x,0) by applications of and, from Corollary 5.1, System 12.1 f |l does not admit an algebraic 
first integral. 

In general, algebraic integrability restricts the possible local behavior of singularities. Theorem 4.3 
gives two possible local behaviors but in each cases, the singularity can have several different global 
behaviors. Consider an irreducible component of the singular set S'^^-'(f). In the proofs of Theorems 
3.3 and 4.3, we considered two cases: S^^\{) (/_ L^ and S^^\{) C Ic- When ^(^^(f) <;t L^ the sin- 
gularities of S'l^-'(f) are generically confined. In this case, two different types of global behavior can 
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be expected. First, after meeting a finite number singularities, the mapping generically never meets 
other singularities (note that this must also be true in the backward direction, that is for f""'^). We 
will referred to such a behavior as global confinement. Second, the singularities may be confined but 
not globally. Because, for each dynamical system of the form (|2.21|l . there are only finitely many 
irreducible parts to 5(f), the mapping must, at one point, come back to the initial singularity and 
take S\^\t) into itself. We say that such singularities are periodic. We now define formally these two 
types of global behavior and illustrate them. 

Definition. A singularity x* G for a birational mapping x„_|_i — f(x„) is said to be globally 
confined if there exist k',k" e such that x* ^ S'(f'^) for any k > k' and x* ^ S{{^^) for any 
k > k". 

An example of such a behavior is given by System (|2.2|l . For this system, points of the form (a;,0) 
are not sent to points of the form (0, y) through successive applications of f . To prove this, let 
{xmUn) be the nth iterate of (x,0). If n = 4fc, an expansion of f*(a;,0) reveals that {x4k,y4k) — 
[kx + 0{x'^), (-fc + l)a; + 0{x^)), k ^ 1,2, . . .. This fact also implies that x^k+z = -kx + 0{x^) 
and we conclude that 0:4^+3 and X4^k are neither infinite nor vanishing. By contradiction assume that 
there exists n such that {xn,yn) = (0,y). Therefore, either n = 4fc + 1 or n = 4fc + 2. However, 
from Equation l|2.4|l we know that Xn+2 is infinite which is a contradiction. Wc conclude that no iter- 
ate of (x, 0) falls onto a point of the form (0, y) and singularities of the form (0, y) are globally confined. 

Definition. Consider the singular set a birational mapping x„_|_i = f(x„) and one of its irreducible 
component ^(^^(f). If there is a k e N>° such that {''{S^^^f)) = S^^^f), then the elements of S^^^f) 
are said to be periodic singularities of period k. 

An example of such a behavior is given by a particular case of the Gambler system |13II14| 



Xn+l 
Vn+l 



f (^ni yn) 



/ iVn + (1/a - a - - a + 1 \ 

V 1 - i/yn / 



(5.5) 



where a is a nonzero complex constant. The system (|5.5|l admits a singularity of type II at any point 
of the form (x, a). Then the first iterate of a point of this form is given by (1/a — 1, 1 — 1/a). But the 
Jacobians of f ^ and f'^ at (x, a) are nonzero and the singularity is confined. It is not globally confined 
since f''(x, a) = f(x, a) = (1/a — 1, 1 — 1/a) and the singularity is periodic. The system admit the 
following first integral 

2/^ - 3y + 1 

I{x,y) = — — . (5.6) 

y{y - 1) 

Now, consider the second case considered in the proofs of Theorems 3.3 and 4.3: 5^^-'(f) C for 
some c. If the singularities of S'^^\t) are generically confined in k steps, then f'^(<S'^^-' (f)) has to be 
a set of the form of the RHS of 1)3. 8|l . Since there are only finitely many of those, the singularities 
cannot be globally confined. So, either they are periodic, or ^^"'^•'(f) lies inside the singular set of f'^ for 
all k greater than a certain k' . It must also be true in the backward direction and if the singularities 
are not periodic S^^\{) lies inside the singular set of for all k greater than a certain k" . We will 
refer to these singularities as ubiquitous singularities. An example of such a behavior is given by the 
singularities of Equation H2.18|l of the form (x, 0). These singularities are not confined but they are 
in the singular set of f ~' for any I > and therefore, these singularities are ubiquitous. Note that 
singularities of the form (0, y) are also ubiquitous. 

Definition. A singularity x* G of a birational mapping is said to be ubiquitous if there ex- 
ist a k' e N>° such that x* G S{{'') for any k > k' and a k" G N>° such that x* G S{i^'') for any 
k > k". 
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Given a discrete dynamical system of the form (|2.21(l and an irreducible component S'(^^(f) of the 
singular set, one would like to be able to obtain information on the first integrals. The only way one 
can know if a given irreducible component of the singular set lies inside a level set is by studying the 
global behavior of the singularity. If the singularity is globally confined, then S''-^''(f) does not lie inside 
a level set, but if it is ubiquitous, then it does. If the singularities are periodic, one cannot reach a 
conclusion. 

Information on first integrals can be obtained by studying singularities of type I. Consider the 
paths Pi defined in the 2-dimensional case in the proof of Theorem 3.1. The number of these paths is 
bounded above by the maximum of the degrees of P and Q. Moreover, for each s'f \{) ^ Lc for all c, 
consider kj, the lowest value kj can take. Then k'^ — 1 is the number of different paths Pi realized by 
iterating a point close to the singularity under the mapping. Hence, the following corollary holds. 

Corollary 5.2. Consider a 2-dimensional hirational mapping with rational first integral R = P/Q. 
Consider also the irreducible components of the singular set of the first type s\^^ (f ) for which the sin- 
gularities are globally confined. Then the maximum of the degrees of P and Q is bounded below by the 
number 

^(fc;.-i), (5.7) 

where k'j is given by the lowest value kj of Theorem 3.1. 

In Example (|2.2|) . the singular set has a unique component in which singularities are globally con- 
fined with fc = 4 and the degree of the first integral is exactly fc — 1 = 3. Note, however, that the 
corollary only provides a lower bound as shown in the following example 

V Vn+l J 

where a is a non-vanishing complex number. Here again, one can check that singularities of the form 
(0, y) are confined in 4 steps and since confinement is global, the lower bound for the numerator and 
denominator of a first integral is 3. However, the rational first integral I = P/Q with lowest degree 
possible on P and Q is given by 

I ^ {Xn~ynY{Xnyn- of, (5.9) 

which is of degree 6, not 3 (note however the rather degenerate nature of the first integral). 

Next, we turn our attention to singularities of type II in the case when S'j{^(f) <t Lc for aU c. 
From the proof of Theorems 3.2 and 4.2, one sees that, as long as the singularities arc not confined, 
the iterates of 5'|{'(f) {Z! Lc must be roots of both Pi and Qi for all rational first integrals Pi/Qi. 
Therefore, we have the following corollary. 

Corollary 5.3. Consider an algebraically integrable p-dimensional birational mapping with the 
rational first integral R = P/Q. Consider also an irreducible component of the singular set of the 
second type Sj''j{f) whose singularities are globally confined. Let kj be the corresponding confinement 
number. Then the finite elements of the set 

kj-l 

U i\S^p{i)). 

are roots of both P and Q. 

This corollary is best illustrated on Example H2.15|l where (0, a) and (1/a, 0) are roots of both the 
numerator and denominator of the rational function (|2.17|) . 
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The previous corollary has an immediate consequence for the existence of polynomial first integrals. 

Corollary 5.4. Consider an algebraically integrable p- dimensional birational mapping. If this system 
admits singularities of the second type which are globally confined, then none of the first integrals are 
polynomial. 

The last result concerns the case when the irreducible component of the singular set lies inside the 
closure of a level set of the first integral. In Sections 3 and 4, it was found that if f'' (for /c S N) is 
well-defined on ^(^'(f) and its Jacobian is nonzero, then f*''(S'(-') (f)) must be a set of the form of the 
RHS of (|3.8|) . hence it lies on a level set. This is stated in a following corollary. 

Corollary 5.5. Consider a p-dimensional birational mapping x„+i — f(x„) together with an ir- 
reducible component of the singular set S^^^{i) for which the singularities are ubiquitous. If f'^ is 
well-defined almost everywhere on S^^\{) and its Jacobian is nonzero, then f''(S^-'\f)) lies inside the 
closure of a level set of the first integrals. Note that the case k = should also be included meaning 
that S^^^{i) itself lies inside the closure of a level set of the first integrals. 

This corollary can be applied to Example (|2.18|l . Since the two sets of singularities {{x,0) and {0,y)) 
are ubiquitous , they are level sets of the first integral (|2.19|l . 
As a final example, consider the system 



where a is a nonzero complex number. This example was first used in to show that singularity 
confinement is not sufficient for integrability. Following Example (|2.2|l . it is easy to show that 
singularities of the form (0, y) are generically globally confined. However, numerical analysis performed 
on this system strongly suggest that is exhibits chaotic behavior [iBj . 

VI Conclusions 

In this paper, we defined different types of singularities and their confinement property for autonomous 
discrete birational dynamical systems. The theory of singularity confinement was developed more than 
10 years ago |llj|(22jll2j and different types of singularities have been implicitly used. Here, in order 
to relate algebraic integrability to singularity confinement we have found necessary to give a formal 
definition of different singularities and their confinement property. 

In Sections III and IV, we showed that algebraic integrability for birational mappings imply that 
singularities are either locally confined or some preimages are not well-defined. 

In Section V, we studied global properties of birational mappings and showed that there are three 
types of behavior for singularities compatible with algebraic integrability. Namely, singularities are 
either globally confined, periodic, or ubiquitous. Remarkably, in the seminal papers describing sin- 
gularity confinement, similar behaviors were found heuristically to be compatible with integrabil- 
ity 13, J^. However, the classical approach relies mostly on a local analysis of the 
singularities. An important aspect of the results shown in this paper is that global properties have to 
be satisfied for a system to be algebraically integrable. To unify these different concepts, we define the 
confinement property for a discrete system in the following way: 

Definition. A birational mapping is said to have the confinement property if its singularities 
are, generically, globally confined, periodic, or ubiquitous. 

Then the most important result of this paper can be rephrased as follows. 

Theorem. An algebraically integrable birational mapping has the confinement property. 
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In particular, this result shows that every mapping in the QRT-family |22II15| has the confinement 
property both in the symmetric and asymmetric cases. This fact also holds for the recently discovered 
class of integrable birational discrete systems discovered in ^18^. Moreover, from examples such as 
System we have strong numerical evidence that the confinement property is necessary but not 

sufficient for algebraic integrability. 

In this paper, we also found sufficient conditions for non-integrability. Essentially, if singularities 
are confined in one direction and not in the other, the system cannot be algebraically integrable. 

Finally, the remaining corollaries of Section 5 provide information on the existence and degree of 
rational first integrals directly from the singularity confinement procedure. 

All the examples of integrable discrete systems considered in this article have a first integral of 
genus or 1. It is not possible to find a birational mapping of infinite order preserving an algebraic 
curve of genus 2 or higher. The argument proving this fact is due to Veselov ( j29j . page 35). It makes 
use of Hurwitz theorem which tells us that the automorphic group of an algebraic curve of genus 2 or 
higher is finite. 

Since the only mappings studied here were birational, we could have made use of the projective 
space in which the point at infinity plays no particular role. Hence only singularities of type II should 
be studied. However, the ideas presented in this paper could not be extended for the case of non- 
autonomous or for discrete dynamical systems that are not birational. For instance, the integrable 
non rational mapping derived in |8IHI1()| could not be studied. It is therefore crucial to give a definition 
of singularity confinement that can be applied to general discrete systems. 

Another important question concerns the relation between the confinement property and other 
types of integrability for discrete dynamical systems (see for example pil51El [TBll24 | ) . 
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